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Abstract 

A semiclassical method of deriving Regge trajectories for the bound 
states of the Klein-Gordon equation with the interaction introduced in 
a non-minimal way is proposed. The method is applied to construction 
of the quarkonium Regge trajectories. It is found that under the rela- 
tivistic generalization of the Cornell potential the Regge trajectories of 
charmonium are in the same good agreement with experimental data for 
introducing the confinement part of potential either in the minimal, or 
non-minimal way. 



1 Intoduction 

Application of potential models is one of the most popular approaches to the 
description of mesons as the bound states of quark- antiquarks systems pQ. Ac- 
cording to QCD predictions, at small distances between the quarks in these 
systems, the Coulomb potential, originating from the one-gluon exchange, dom- 
inates. With increasing the distance, the forces arise which provide the quark 
confinement. In approximation of one-«dressed»-gluon exchange, at large dis- 
tances we have the linear behaviour of the potential that is a theoretical ground 
for introducing the Cornell potential, V(r) = —a/r + nr, in the non-relativistic 
approach PJ. 

Under the generalization to the relativistic case, within the framework of the 
model of quasi-independent quarks, one usually considers the Dirac or Klein- 
Gordon equations. Then, to avoid the Klein paradox, the linear part of the 
potential in these equations is traditionally included in the Lorentz-scalar com- 
ponent PJE]- However, a «non-minimal> inclusion of the linear interaction 
attracts sufficient attention nowadays. The Dirac equation with such an inter- 
action was independently considered by different authors [SH^IE]. But the main 
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investigations were stimulated by the work |Hj after which this model is referred 
to as the Dirac oscillator |§] . For this model, the symmetric and supersym- 
metric properties are studied, the one-dimensional and two-dimensional 
cases are considered. The oscillator interaction is also generalized to the rel- 
ativistic Klein-Gordon ^1 ^5] and Dufhn-Kemmer-Petiau ^] equations. 
However, all above investigations deal with the only unperturbed oscillator in- 
teraction. 

The goal of the present work is to develop the method of deriving Regge 
trajectories for the bound states of the Klein-Gordon equation with the interac- 
tion introduced in the non-minimal way supplemented with the Lorentz-scalar 
potential and the time component of the Lorentz vector. The obtained Regge 
trajectories are applied to study the properties of charmonium mass spectrum 
and to find the difference in the relativistic generalization of the Cornell potential 
with inclusion of the confinement interaction in the minimal and non-minimal 
ways. 



2 Recurrent relations for Regge trajectories 

A general receipt of non-minimal inclusion of the oscillator interaction consists 
in the replacement of the momentum, in the Hamiltonian of a free particle, by 
the rule p — > p — i(3W{r) where f3 is a matrix which in the case of the Dirac 
equation equals to {3 = 70. 

We will study the properties of the non-minimal interaction for a spinless 
particle. This restriction is justified by the fact that the spin effects are not 
essential in quarkonia spectra which we intend to describe. 

Introducing the non-minimal interaction in the Hamiltonian of a spinless 
particle, written in the Sakata-Taketani form, with W(r) = kt yields the radial 
equation for the Klein-Gordon oscillator |121 113] 
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This equation has the exact solution and infinitely degenerate energy spectrum 

E 2 — m 2 — Kh(4n + 21), n, I = 0, 1, 2, . . . (2) 

where n and I are the radial and orbital quantum numbers, respectively, and we 
put c = 1 . 

For eliminating the unphysical degeneracy, the perturbation should be added 
to the oscillator interaction. Let us consider the general case of the radial Klein- 
Gordon equation for a particle of mass m moving in the field of central potential 
consisted of the non-minimal part, W(r) = W(r)f/r, the Lorentz-scalar part, 
5(r), and the time component of the Lorentz vector, V(r): 



tfu"{r) = [ % -{E- V(r)f + W 2 (r) 



-n(W'(r) + -W(r)) + (m + S(r)) 2 ]u(r). (3) 
r 
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The potential functions S(r), V(r), W{r) are assumed to satisfy the condition 
that the effective potential has two turning points and the equation leads to 
discrete energy spectrum. 

Our task is to derive the Regge trajectories, a(E) = M(E), for the bound 
states of Eq. © . For this end, we generalize the semiclassical technique proposed 
in the works |16II17| to the case of the non-minimal interaction. 

With the substitution C(r) = hu'(r)/u(r), we rewrite Eq.© in the form of 
the non-linear Riccati equation 

hC'(r) + C 2 (r) = ^a{E){a(E) + h) — (E — V(r)f 

+(m + S(r)) 2 + W 2 (r) - h(W'(r) + -W(r)). (4) 

r 

We will seek its solution as the series expansion in powers of Planck's constant 

oo oo 

C(r) =Y / C k (r)h k , a{E) = M{E) = Y d <Xk{E)h k . ( 5 ) 

k=0 fc=0 

Substituting the expansions (J5J into Eq. and collecting the coefficients at the 
same powers in h, we have 

C 2 (r) = a 2 {E)/r 2 - (E - V{r)f + (to + S{r)) 2 + W 2 (r), 

j=0 3=0 

-S kA (W'(r) + -W(r)), k>l. (6) 
r 

In order to take into account the radial quantum number n in the quantum 
corrections, we apply the argument principle, known from the complex analysis, 
to the logarithmic derivative C(r) that results in the following quantization 
condition 

-±-.j>C{r)dr = nh, 71 = 0,1,2,..., (7) 

where n is the number of nodes of the wave function on the real axis, and the 
integration contour encloses only the above nodes. 

This quantization condition should be supplemented with a rule of passing 
to the classical limit for the radial and orbital quantum numbers which are the 
specific quantum notions. Recall that in the Wentzel-Kramers-Brillouin method 
the passage to the classical limit is realized by the rule 

h — ► 0, n — > oo, I — > oo, hn = const, Til = const. (8) 

For deriving the Regge trajectories, we use an alternative opportunity 

h — ► 0, n = const, I — > oo, hn — > 0, M — const (9) 
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that permits us to rewrite Eq.lQ) in the form 



■^l <f Ci(r)dr = n, ^jc k {r)dr =0, k ± 1. (10) 

From the physical point of view, such passage to the classical limit corresponds 
to the situation when the particle is located at the bottom of the potential well 
and thus moves along the stable circular orbit of radius r . The particle has the 
angular momentum ao and the energy 

E=^J{m + S(r Q )) 2 + W 2 (r Q ) + a 2 /r 2 + V(r ). (11) 

Then it must be Co(tq) — 0, and the radius of orbit is calculated as the position 
of the minimum of effective potential from the equation 

ro[V'(r )(E - V(r )) + S'(r )(m + S(r Q )) + W'(r )W(r Q )} 

= [E- V(r Q )) 2 - (m + 5(r )) 2 - W 2 (r ). (12) 

Consequently, for the classical angular momentum ao(E) which is the first 
approximation to the Regge trajectory, we have 

a (E) = ^Jr 3 [V'(r )(E - V(r )) + S'(r )(m + S(r )) + W'(r )W(r Q )}. (13) 

Further, it is convenient to introduce a new variable, x = (r — ro)/ro, that is 
the deviation from the minimum of effective potential. Because Co(ro) = 0, the 
function Cq(x) can be represented in the form 

oo 

Cq(x) = -UJXy^l + aix + a 2 x 2 + ... = x^C°x\ (14) 

i=0 

where C° = C? = i(£ tfC^ ~ u 2 ^), i > 1. 

3=1 

Coefficients u>, at, in their turn, are expressed through the coefficients of 
expansions of the potentials in the Taylor series in powers of x 

oo oo oo 

V(r) = S(r) = S t x\ W(r) = £ W t x l (15) 



as 



uj 2 = ial/rl + 2{E - V )V 2 - V? + 2(m + S )S 2 + S 2 + 2W W 2 + W?, 
a, = [(-l) l (3 + i)a 2 /r 2 + 2EV i+2 + 2mS l+2 

i+2 

+ /J (SpSi+ 2 - p — V p Vi+ 2 -p + WpW. l+2 -p))]/u 2 . (16) 



p=0 
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Consider now the functions Ck(x), k > 1. Since Co(0) = 0, from the struc- 
ture of Eqs.JBJ) it follows that the function Ck(x) has the pole of (2k — l)-th 
order in the point x = and can be expanded into the Laurent series 



oc 



C k (x)=x 1 - 2k J2^x i . (17) 



i=0 



Such representation of the functions enables, by using the residue theorem, 
express the quantization conditions (|10|l directly through the coefficients C\ 

C\ = n/ro, C* k _ 2 = 0, k ? 1. (18) 

Eas. (|14|) and l|18|l permit us to compute the part of coefficients Of. Substituting 
the expansions l|14|l and l|17|) in Eqs.|jf)| yields the expression for C* with k > 1 
and i ^ 2k - 2 

1 , i 2k ■ 3 ''" L ' ' 



k 

3 = 1 P=0 P=l 

(-l) 4 (i-2fc + 3) 



+0(i - 2k + 2)P — g -(a*-i + ^ Ojafc-j-) + E Oj-afc-j-) 

1 ' 
-«k,i — ((* + + 2 ^ (-1)%-,)]}. (19) 

r ° p=0 

Whereas in the case of i — 2k — 2 we arrive at the recurrent relation for the 
coefficients of ^-expansions for the desired Regge trajectories 

1 fc_1 

a k (E) = --—[a k -i + E a J a k-j - r o c %k-i 
2a 

k 2k-2 

~ r o E E C l C 2k-2- P - h,MWi + 2W )]. (20) 

j=0 p=0 

The obtained recurrent formulae completely resolve the problem of deriving 
the Regge trajectories for the bound states of the radial Klein-Gordon equation 
with the non-minimal confinement interaction supplemented with the Lorentz- 
scalar and Lorentz- vector potentials. 



3 Charmonium Regge trajectories 

In order to check the possibility of the non-minimal inclusion of the quark 
confinement interaction, let us apply the developed approach to the description 
of the charmonium mass spectrum. 

It should be stressed, that the relativistic generalization of the Cornell po- 
tential has enough arbitrariness. For example, in the string model the account of 
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quantum fluctuations leads to appearence of the Coulomb addition to the long- 
range potential ^3 EH which is assumed to be the Lorentz-scalar. Whereas 
in the non-relativistic problem this addition is combined with the short-range 
Coulomb potential, now we are faced with the problem of redistribution of the 
Coulomb part between Lorentz- vector and Lorentz-scalar components. There- 
fore in our model we choose the potentials in the form 

W(r) = nr, V(r) = -b/r, S(r) = -a/r. (21) 

For these potentials, restricting to the second order in h in the expansions 
for Regge trajectories that provides a sufficient accuracy for practical purposes, 
from Eq. l|20|) we obtain 

a( E ) = «r yr+7+^+M *"^^ -h + 

2vl + p + (J 2 

-2 



rr r (3-<h/4 + pr 



2kt^1 + p + a 4(1 + p + cr) 

{2p 2 + 25p + 128)g 2 + 65p + 288 q(3p + 48) 
+ 8(p + 4) 2 2(p + 4) 3 / 2 J 



(22) 



where p = (Eb + ma)/ K 2 r$, a = (b 2 — a 2 )/K 2 r^ q = 2n + l and r$ is a positive 
root of the equation 2k 2 Tq — (E 2 — rn 2 )r — (bE + am) = 0. 

To adopt the formula 122(1 to the description of charmonium, we must sub- 
stitute the quark masses, m — m c = m £ , and express the energy through the 
charmonium mass as E = M/2. 

The minimization of existing experimental data |20| results in the following 
values of parameters (in units with h = c = 1) 

7ti c = 2M2GeV, k = 0.1210(GeF) 2 , a = 1.687, b = 0.051. (23) 

For comparison, we also calculated the Regge trajectories in the conventional 
potential model with 

W(r) = 0, V(r) = -b/r, S(r) = nr - a/r. (24) 

In this case, the best description is achieved at the values 

m c = 1.895GeV, k = 0.0424(GeF) 2 , a = 1.666, b = 0.093. (25) 

Table 1 presents the charmonium mass spectrum calculated within the frame- 
work of these two models. 

Behaviour of the Regge trajectories is shown in Fig.l where solid lines cor- 
responds to the model (|21|) with the non-minimal interaction and dotted ones 
to the conventional approach 124|) . 

From Fig.l we see that the inclusion of the confinement part of the potential 
in the scalar component or in the non-minimal way does not principally change 
the behaviour of charmonium Regge trajectories. 
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Table 1: Charmonium mass spectrum. 



State 


Mass (MeV) 


Experiment 


Model J2JJ with non- 
minimal interaction 


Model 424ft 


J/V 


IS 


3096.87±0.004 


3096.90 


3096.86 




IP 


3556.18±0.13 


3556.35 


3556.59 


V>(3685) 


2S 


3685.96±0.09 


3685.62 


3685.56 


V>(4040) 


3S 


4040±10 


4041.85 


4043.61 



Figure 1: Charmonium Regge trajectories 
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It has also been found that the Regge trajectories in the both models un- 
dergo qualitatively similiar changes under the redistribution of Coulomb part 
of the potential between the Lorentz-scalar and Lorentz-vector. Note that the 
inclusion of extra Coulomb potential in the non-minimal way leads to worse 
description. Indeed, in this case the Cornell potential is not restored in the 
non-relativistic limit. 

Thus, we may conclude that Regge trajectories of charmonium are in the 
same good agreement with experimental data for introducing the confinement 
part of potential either in the minimal, or non-minimal way. The final decision 
in this matter requires additional experimental data on the orbitally-excited 
charmonium states, especially, for the states lying on the daughter trajectories. 
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